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Abstract. We have used a detailed non-LTE radiative transfer code to model new APEX CO(7 = 3 — » 2) data, and existing CO 
radio line data, on a sample of 40 AGB S-stars. The derived mass-loss-rate distribution has a median value of 2 x 10~ 7 M Q yr _1 , 
and resembles values obtained for similar samples of M-stars and carbon stars. Possibly, there is a scarcity of high-mass-loss- 
rate (> 10~ 5 M Q yr~') S-stars. The distribution of envelope gas expansion velocities is similar to that of the M-stars, the median 
is 7.5 kms~', while the carbon stars, in general, have higher gas expansion velocities. The mass-loss rate correlates well with 
the gas expansion velocity, in accordance with results for M-stars and carbon stars. 
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1. Introduction 

Stars of spectral class S are identified through their strong 
bands of primarily ZrO. This indicates enhancements of s- 
process elements. It is now relatively well-established that 
these enhancements can be produced in two ways: through 
dredge-up of nuclear-processed material during the asymp- 
totic giant branch (AGB) evolution (intrinsic S-star), or through 
mass transfer in a binary system (extrinsic S-star). The strong 
bands of ZrO in intrinsic S-stars, as opposed to the strong bands 
of TiO in M-stars, have been attributed to a stellar atmosphere 
chemistry in a relatively low-temperature (^3000K) gas en- 
riche d in s-process elements and with C/O = 1 (w ithin about 
5%) JScalo A ROSSII1974 [Smith A T,arnherll ll990l). This then 
lead to the belief that the AGB S-stars are transition objects 
between AGB M-stars and carbon stars. 

Extensive mass loss is a key process during the evolution of 
AGB stars. M-stars and carbon stars lose substantial amounts 
of matter during this phase, and th eir mass-loss cha racteristics 
are relatively well-established (see Olofssonl l2003h . The puta- 
tive transition objects, the S-stars, have attracted less attention 
in this respect. Studies of CO radio line and dust continuum 
emission suggest that they behave in roughly the same way as 
both their predecessors and descendants in terms of mass-loss 
rate, circ umstellar kinemati cs, dust co mpositio n, and dust-to- 
gas ratio J.Turalll988t iBieging & Latteill 19941 ISahai & Liechtil 
1995l lGroenewegen & de .Tongll998H.Torissen & Knappll9 98). 



In this Letter new APEX 1 CO(7 = 3 -> 2) data, in combi- 
nation with existing CO radio line data and a radiative transfer 
code, were used to derive the circumstellar properties of 40 S- 
stars on the AGB. 



2. Observations 

The sample was selected from the list of 65 objects with S -star 
charac teristics in the Two Mic ron Sky Survey bv lWing & Yorkal 
dl977l) [see also [jura] Jl988l) l. Observations of CO J = 3—>2 
line emission at 345.80GHz towards a sample of 18 S-stars 
were performed from August to October 2005 using the APEX 
12 m telescope equipped with a double-sideband SIS receiver 
and a FFTS spectrometer. The system temperature (above at- 
mosphere) varied between 130 and 335 K. The effective band- 
width at the observations were 1 GHz with a spectral resolu- 
tion of 61 kHz (0.05kms _1 ), except for FU Mon, R Gem, and 
Z Ant observed during the second run with a spectral resolution 
of 122kHz (O.llkms -1 ). The data was smoothed to lkms -1 
resolution and reduced by removing a low-order (< 3) polyno- 
mial baseline in the range of about + 100 km/s from the stellar 
velocity. The resulting spectra are presented in Figure. ^ an d 
velocity-integrated intensities are reported in Tabled 

The observations were carried out using a position- 
switching mode, with the reference position located +2' in az- 
imuth. The raw spectra are stored in the T* scale and converted 
to main-beam brightness temperature using T m \, — T^/j] m b, 
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1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Fig. 1. New observations of CO(/ = 3— > 2) line emission using the APEX telescope. The velocity resolution is l.Okms 



where T* is the antenna temperature corrected for atmospheric 
attenuation using the chopper-wheel method, and rj m \, the main- 
beam efficiency, estimated to be 0.7 at 345 GHz. The uncer- 
tainty of the absolute intensity scale is estimated to be about 
+ 20%. Regular pointing checks were made on strong CO 
sources and typically found to be consistent within ss 3". The 
FWHM of the main beam, mb , is 18" at 345 GHz. 

In addition to the new APEX data presented here, 
we have also used new CO(7 = 1 — > 0) data observed 
at the Onsala 20 m telescope ( January 2006), JCMT 2 
public ar chive CO data, and C O line intensi ties re- 
ported bv iBieeing & LatteJ Jl994. ISahai & Liechtil jl 995). 
Staneket alJ lll995h. iGroenewegen & de Jong l Jl998l) . and 
d 1JJ998). 



Jorissen & Knap 



In total, 42 S-stars have been de- 



tected in CO line emission. 



3. Radiative transfer analysis 

In our analysis the circumstellar envelopes (CSEs) are assumed 
to be spherically symmetric, to be produced by a constant mass- 

2 The James Clerk Maxwell Telescope, is operated by the Joint 
Astronomy Centre in Hilo, Hawaii on behalf of the parent organiza- 
tions PPARC in the United Kingdom, the National Research Council 
of Canada and The Netherlands Organization for Scientific Research. 



loss rate (M), and to expand at a constant velocity (v e ). The 
spectra of n l Gru and RS Cnc indicate that deviations from 
the standard CSE model are present and no mod e lling was at- 
tempted for these two sources (see ISahailll992t iKnapp et alJ 
1999). 



3.1. SED modelling 

The SEDs were constructed from IRAS-fluxes and 2MASS- 
magnitudes. The continuum emission was modelled us ing th e 
dust radiative transfer code DUSTY (Ivezic & Elitzurlfl99l) . 
Amorphous carbon or amorphous silicate dust was selected 
based on, in most cases, the sp ectral class of the IR AS low- 
resolution spectra as defined bv IVolk & Cohenl dl989l) . In the 
cases where the spectral class was not available, the selec- 
tion was based on col our-colour diagrams as described by 
IJorissen'&' Knapp (1998). The adopted dust type is listed in 
Tabled where C and O denotes carbon and silicate dust, re- 
spectively. The amount of dust present in the wind can be con- 
strained from the spectral shape of t he SED. Details on the 
modellin g procedure c an be found in Sch oier et alJ (|2002j) and 
Ram stedt et all feOOd l. 

The SED results were used in combination with the lu- 
minosities to estimate the distances to the sample stars. 
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Table 1. APEX observations and model summary 



CO modelling 



Source 


D 


u 


7(3 - 2) 


M 




h 


*r 2 ed 


N 


dust 




[pc] 


[L ] 


[Kkms>] 


[M s yr-'] 


[kms- 1 ] 








type 


R And 4 


300 


6000 




6.0 x 10~ 7 


8.0 


0.2 


0.1 


2 


O 


W And" 


280 


5800 




1.2 x 10~ 7 


5.5 


0.2 


1.0 


3 


O 


Z Ant 


470 


4000 


2.0 


9.0 x 10~ 8 


6.0 


0.2 




1 


O 


VX Aql 


790 


4000 




5.0 x 10~ 8 


7.0: 


0.2 




1 


C 


W Aql c 


230 


6800 


136 


2.5 x 10~ 6 


17.5 


0.7 


2.5 


4 


O 


AA Cam 


800 


4000 




1.8 x 10~ 7 


16.9: 


0.2 




1 


C 


T Cam 


540 


5600 




1.0 x 10~ 7 


4.8: 


0.2 




1 


C 


S Cas 


440 


8000 




2.2 x 10~ 6 


19.0 


0.5 


0.9 


2 


O 


V365 Cas 


625 


4000 




3.0 x 10~ 8 


6.2: 


0.2 




1 


C 


WY Cas 


600 


6700 




6.0 x 10 7 


12.5: 


0.5 




1 


O 


AM Cen 


750 


4000 


1.1 


1.5 x 10~ 7 


3.3 


0.2 


1.9 


2 


C 


TTCen 


880 


6500 


13.9 


2.5 x 10~ 6 


20.0 


0.5 


7.8 


2 


O 


UYCen 


590 


4000 


1.6 


1.2 x 10~ 7 


12.0 


0.2 




1 


C 


V386 Cep 


470 


4000 




1.8 x 10~ 7 


18.6: 


0.2 




1 


O 


T Csi aicd 


240 


4000 


3.7 


5.0 x 10~ 8 


5.5 


0.2 


0.8 


4 


C 


AA Cyg 


480 


4000 




4.0 x 10~ 7 


6.1: 


0.2 




1 


C 


RCyg" 


440 


6100 




5.0 x 10~ 7 


9.5 


0.5 


1.8 


3 


O 




1 10 


5900 




5.0 x 10~ 7 


8.5 


0.1 


2.0 


3 


O 


TV Dra" 


390 


4000 




7.0 x 10~ 8 


4.0 


0.2 




1 


O 


DY Gem 


680 


4000 




1.3 x 10~ 7 


7.0 


0.2 




1 


O 


R Gem" 


710 


5500 


6.4 


4.0 x 10~ 7 


4.0 


0.5 


1.0 


2 


C 


ST Her" 


300 


4000 




1.3 x 10~ 7 


8.5 


0.2 


1.1 


3 


O 


RX Lac 


310 


4000 




6.5 x 10~ 8 


5.5: 


0.2 




1 


C 


GILup 


690 


5000 


4.6 


1.0 x 10~ 6 


10.0 


0.2 




1 


C 


RLyn" 


850 


5600 




3.0 x 10~ 7 


7.0 


0.2 


0.2 


2 


C 


YLyn" 


260 


4000 




1.7 x 10~ 7 


7.0 


0.2 




1 


O 


S Lyr 


1210 


6300 




1.0 x 10~ 6 


13.0: 


0.5 




1 


O 


FU Mon 


780 


4000 


1.3 


1.2 x 10~ 7 


2.0 


0.2 




1 


c 


RZ Peg 


970 


4000 




3.4 x 10~ 7 


12.6: 


0.2 




1 


c 


RTSco 


270 


6400 


32.6 


4.5 x 10~ 7 


11.0 


0.5 


3.2 


2 


o 


ST Sco 


380 


4000 


4.5 


1.3 x 10~ 7 


5.5 


0.2 


1.2 


3 


o 


RZSgr 


730 


4000 


15.8 


1.5 x 10~ 5 


9.0 


0.2 


22.3 


2 


c 


ST Sgr 


540 


5800 


4.9 


2.7 x 10~ 7 


7.0 


0.2 


3.0 


2 


o 


TSgr 


590 


5800 


2.5 


2.0 x 10~ 7 


10.0 


0.2 




1 


c 


EPVul 


510 


4000 




1.5 x 10~ 7 


4.4: 


0.2 


5.5 


2 


c 


DK Vul 


750 


4000 




6.0 x 10~ 7 


4.0: 


0.2 




1 


c 


AFGL2425 


610 


4000 


2.8 


2.0 x 10~ 7 


6.5 


0.2 


0.9 


2 


o 


CSS2 41 


880 


4000 




6.0 x 10- 7 


20.0: 


0.2 




1 


o 


IRC-30285 


590 


4000 


6.5 


6.0 x 10~ 7 


17.0 


0.2 


6.0 


2 


o 


IRC- 10401 


430 


4000 


9.5 


3.5 x 10~ 7 


17.0 


0.2 


12.8 


2 


o 



Superscripts mark new or archive data: " OSO y = 1 — > 0, 
b JCMT J = 3 -> 2 , c JCMT J = 4 -> 3, d JCMT 7 = 6^5. 
The forth column gives the APEX J = 3 — > 2 results. 



Luminosities of the M iras were estim a ted us ing the period- 
luminosity relation of Whit elock et alJ Jl994l) . For the semi- 
regular and irregular variabl es, a luminosity of 40 00 Lp was as- 
sumed in accordance with lOlofsson et alJ ( 120021) . Corrections 
for interstellar extinction were applied. The derived distances 
and adopted luminosities are listed in Tabled Non-negative 
Hipparcos parallaxes with an accuracy better than 50 % are 
available for seven of the sample stars. They agree well within 
the error margins, when compared to the distances derived from 
the SED modelling (except for RZ Peg). 

The results from the SED modelling were also used to de- 
termine the dust radiation field which is important for the exci- 
tation of CO. 



3.2. Radiative transfer model 

In order to determine the molecular excitation in the CSE, a 
detailed non-LTE radiative transfer code, based on the Monte 
Carlo method, was used. This code is described in detail in 
Schoier & Olofsson <200 111 . It has been benchmarked, to high 
accuracy, again st a wide variety of molecular-line radiative 
transfer codes in van Zadelhof f et alJll2002l) . 

The CO molecules are partly excited through c ollisions 
with H ?. The collisional rate coefficients calculated by Flower 
fcOOlh were adopted and extended to include more energy lev- 
elsas well as extrapolated in temperature (see ISchoier et alJ 
2005). An ortho-to-para ratio of 3 was adopted when weight- 
ing together collisional rate coefficients for CO in collisions 
with ortho- and para-H2. The excitation analysis includes ra- 
diative excitation through the first vibrationally excited (v= 1) 
state at 4.6 jum. Rel evant molecular data are summarised in 
ISchoier et alJ d2005t) and are made publicly available through 
the Leiden Atomic and Molecular Database (LAMDA) 3 . A 
central blackbody was assumed to represent the stellar radi- 
ation field, and thermal, circumstellar dust emission was in- 
cluded. These provide the main sources of the infrared photons 
that excite the v = 1 state. The addition of a dust component 

in the Monte Carlo scheme is straightforward as described in 

■ 1 j — ■ — -j a 

Schoier et al. (2002). The energy-balance equation was solved 
simultaneously as the excitation of CO, allowing for an accu- 
rate determination of the kinetic gas temperature of the gas. The 
photospheric abundance of CO, relative to H2, was assumed to 
be 6x 10~ 4 . A micro-turbulent velocity of 1 km s _I was as- 
sumed. 

Two parameters are adjustable when fitting the model to 
the observed line intensities: the mass-loss rate M and the h- 
paramete r. The latte r is a product of different dust parameters 
( Sch oier et all2002l) . Collisional heating through gas colliding 
with dust grains is the dominant heating term in the energy- 
balance equation and it affects the model line-intensity ratios. 
We were able to constrain h for five of our stars. Its dependence 
on stellar luminosity is consistent with our findi ngs for M-stars 
(lOlofsson et alJ Eo02l) and carbon stars ( Sch oier & Olofssonl 
2001). Therefore, a value was assumed for the remaining stars 
depending on the stellar luminosity, h-0.2 for L < 5000 L and 
h-0.5 above 5000 L (indicated by a colon in TableQ. 

The best-fit model was found by minimizing the difference 
between the observed and modelled integrated line intensities 
(N, in Tabled is the number of observed lines used). Reduced 
^-values were calculated and are presented in Tabled when a 
sufficient number of lines were available. The uncertainties in 
the derived mass-loss rates were estimated to be about ± 50% 
within the adopted circumstell ar model, excluding distance un- 
certainties, in accorda nce with lSchoier & Olofssonl 12001 ) and 
lOlofsson et all J2002I) . 

A particular problem occurred with the CO(J = 2 — > 1 ) data. 
We found a large scatter in the reported intensities for an in- 
dividual object, even when observed with the same telescope. 
Furthermore, these data are, in general, not compatible with 
the model results, i.e., the 3 — > 2/1 — » intensity ratios are 

3 http : //www. strw. leidenuniv.nl/~moldata 
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Fig. 2. a) Mass-loss-rate distributions for the S-star (solid line; 40 stars), M-star (dotted line; 77 stars), and carbon star (dashed 
line; 61 stars) samples, b) Envelope gas expansion velocity distributions for the S-star (solid line), M-star (dotted line), and 
carbon star (dashed line) samples, c) Derived mass-loss rates plotted against the gas expansion velocities for the S-star (dots), 
M-star (squares), and carbon star (triangles) samples. 



compatible with the model results, while the 2 — > 1/1 — >0 or 
3 — > 2/2 — > 1 intensity ratios are, in general, not (note, the model 
line-intensity ratios only weakly depend on the assumed pa- 
rameters). Consequently, we decided to omit the 2 — > 1 data 
from our analysis, except for the few cases where only a 2 — > 1 
line exists. 



lar model, and accurate molecular abundances will be derived. 
These will be compared with results of chemical modelling, as 
well as with the results for the M-star and carbon star samples. 

Acknowledgements. SR, FLS, and HO acknowledge financial support 
from the Swedish Research Council. 



4. Results and discussion 

We compare here our results for the S-star sample, which 
were obtained through detailed radiative transfer modelling, 
with t hose obtained in the same way for a sample of M- 
stars Joi ofsson et al. 120021) and a sample of carbon stars 
JSchoier & Ok)fesonf 200 lb . The samples are flux-limited but 
not complete, except for the carbon-star sample, which con- 
tains a sub sample that is most likely c omplete out to a distance 
of 500 pc JSchoier & Olofssor]l200lh . The incompleteness of 
the samples and the small number statistics limit the generality 
of our conclusions. 

Figure^ shows the mass-loss-rate distribution of the S-star 
sample in comparison with those of the M-stars and the car- 
bon stars. The S-star mass-loss-rate distribution appears very 
similar to that of the M-stars and the carbon stars. The me- 
dian value is 2 x 1O~ 7 M yr _I for the M- and S-star samples, 
and 3 x 10~ 7 M Q yr~' for the carbon-star sample. There is pos- 
sibly fewer S-stars with a high mass-loss rate (> 10~ 5 M yr -1 ). 
Figure^ shows the expansion velocity distribution of the S- 
star sample in comparison with those of the M-stars and car- 
bon star samples. There is an indication that the carbon stars 
have higher expansion velocities, while those of the S- and M- 
stars appear similar. The median expansion velocity is 7.5 km 
s _1 for the S-star and M-star samples, and 11km s for the 
carbon star sample. Figure^ shows that the relations between 
mass-loss rate and envelope expansion velocity are very simi- 
lar for the three samples, pointing towards the same mass-loss 
mechanism. 

This is the first step in a project aimed at obtaining an ex- 
tensive data base of circumstellar molecular lines towards a 
sample of S-stars. The CO data will constrain the circumstel- 



References 

Bieging, J. H. & Latter, W. B. 1994, ApJ, 422, 765 
Flower, D. R. 2001, J. Phys. B., 34, 2731 
Groenewegen, M. A. T. & de Jong, T. 1998, A&A, 337, 797 
Ivezic, Z. & Elitzur, M. 1997, MNRAS, 287, 799 
Jorissen, A. & Knapp, G. R. 1998, A&AS, 129, 363 
Jura, M. 1988, ApJS, 66, 33 

Knapp, G. R., Young, K., & Crosas, M. 1999, A&A, 346, 175 
Olofsson, H. 2003, in Asymptotic giant branch stars, ed. H. J. 

Habing & H. Olofsson (Astronomy and astrophysics library, 

New York, Berlin: Springer), 325 
Olofsson, H., Gonzalez Delgado, D., Kerschbaum, F, & 

Schoier, F. 2002, A&A, 391, 1053 
Ramstedt, S., Schoier, F. L., Olofsson, H., & Lundgren, A. A. 

2006, A&A, to be submitted 
Sahai, R. 1992, A&A, 253, L33 
Sahai, R. & Liechti, S. 1995, A&A, 293, 198 
Scalo, J. M. & Ross, J. E. 1976, A&A, 48, 219 
Schoier, F. L. & Olofsson, H. 2001, A&A, 368, 969 
Schoier, F. L., Ryde, N., & Olofsson, H. 2002, A&A, 391, 577 
Schoier, F. L., van der Tak, F. F. S., van Dishoeck, E. F, & 

Black, J. H. 2005, A&A, 432, 369 
Smith, V. V. & Lambert, D. L. 1990, ApJS, 72, 387 
Stanek, K. Z., Knapp, G. R., Young, K, & Phillips, T. G. 1995, 

ApJS, 100, 169 

van Zadelhoff, G.-J., Dullemond, C. P., van der Tak, F. F. S., 

et al. 2002, A&A, 395, 373 
Volk, K. & Cohen, M. 1989, AJ, 98, 931 
Whitelock, P., Menzies, J., Feast, M., et al. 1994, MNRAS, 267, 

711 

Wing, R. F. & Yorka, S. B. 1977, MNRAS, 178, 383 



